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ABSTRACT 

We present two epoch polarimetric images of the BL Lac object 0735+178 obtained with the Very Long Base- 
line Array at 22 and 43 GHz. These images provide the highest resolution observations of this source to date, and 
show a twisted jet with two sharp apparent bends of 90° within two milliarcseconds of the core. The magnetic 
field appears to smoothly follow one of the bends in the jet, suggesting that this structure may be the result of a 
precessing nozzle in the jet of 0735+178. Proper motion information is not conclusive as to whether the twisted 
geometry is the result of knots moving ballistically or following the direction of the jet axis. 

Subject headings: Polarization - Techniques: interferometric - galaxies: active - BL Lacertae objects: individual: 
0735+178 - Galaxies: jets - Radio continuum: galaxies 



1. INTRODUCTION 

0735+178 was fi rst classified as a BL Lacertae object by Car- 
swell et al. (1974), who identified a pair of absorption lines in 
, an otherwise featureless optical spectrum, thereby obtaining a 
lower limit to the redshift of the source z m i n = 0.424. Hence, 
any observed apparent speeds in 0735+178 calculated using this 
redshift are strictly lower limits. 

Radio maps of this source at milliarcsecond resolution have 
been obtained by using VLBI arrays at centimeter wavelengths, 
and show a compact core and a jet of emission extending to 
the northeast. Polarimetric VLBI observations of this source 



were first performed by Gabuzda, Wardle & Roberts (1989) at 
a wavelength of 6 cm: the corresponding polarized intensity 
image shows a jet magnetic field that is predominantly perpen- 
dicular to the jet axis, together with a dramatic change in the po- 
larized flux by 40% over a period of 24 hours. The multi-epoch 
6 cm V LBI o bservations of 0735+178 presented by Baath & 
Zhang (1991) indicated superluminal motion in their compo- 
nent CO at an apparent velocity of ~7.9 h~ l c (H = 100 h km 
s _1 Mpc _I , q = 0.5), moving between the core and a station- 
ary component, which they labeled B, situated about 4.2 mas 
from the core. This situa tion re sembles that found in the quasar 
4C 39.25 (Alberdi et al. 1993 and references therein), where a 



very strong component moves superluminally between the core 
and an outer stationary component, and its increasing bright- 
ness as it moves is interpreted as Doppler enhanc ement c aused 
by a bend toward the line of sight. Gabuzda et al. (1994a; here- 



after G94) confirmed the motion of the superluminal compo- 
nent CO, which they designated as K2, at a velocity of 7.4 h~ l c, 
and detected two new superluminal components with apparent 
transverse speeds of 5.0 h~ l c and 4.2 h c. The observations 
of G94 were near the time of intersection of the moving com- 
ponent K2 with the stationary one Kl (component B of Baath 
& Zhang 1991 ), which took place at epoch ^1989.8. Their im- 
ages do not show any evidence for a violent interaction between 



these components, and they interpret this in t erms o f a model 
similar to that suggested by Baath & Zhang (1991), in which 
the stationary component ,K7 is associated with a bend in the jet 
toward the line of sight. However, if this is the case, we might 
expect a deceleration of K2 as it approaches Kl, along with an 
increase of its flux due to enhancement of Doppler boosting, as 
observed and simulated in the case of 4C 39.25 (Alberdi et al. 



1993). Neither of these were observed in 0735+178. If Kl cor- 



responds to a bend toward the observer, it also remains unclear 
why G94 measured component Kl to have shifted ~ 0.7 mas 
outward sometime between 1981.9 and 1983.5, after its colli- 
sion with K2, and to have remained in that position afterward. 
One explanation for this "dragging" of Kl is that we are dealing 
with an interaction between a moving sh ock (K 2) and a stand- 
ing conical shock (Kl), as Gomez et al. (1997) have modeled 
through numerical simulations. In this scenario, an increase in 
the Mach number of the flow causes the stationary shock to re- 
set to a new site downstream of its original position. 



Higher resolution VLBI maps obtained by Baath et al. ( 199 1 ) 
at 1.3 cm reveal a complex structure consisting of several dis- 
tinct superluminal components moving in a rather straight jet 
extending to the northeast. However, a fit to the kinemati- 
cal pro perties of the superluminal components led Baath et al. 
(1991) to suggest a bend of the jet in the inner milliarcsecond. 
Indeed, a quite complex structure in the inner milliarcsecond 
of 0735+178 ap pears in the 1.3 cm VLBI maps presented by 
Zhang & Baath (1991). Over a time span of about 1.2 yrs, the 
jet structural position angle was found to chang e from ~ 25° to 
~ 73°. It is unclear whether Zhang & Baath (1991) observed 
different components tracing a bent common path, or compo- 
nents being ejected from the core with different position angles, 
as G94 suggested based on their observations. The first direct 
evidence of a curved structure in the in ner jet of 0735+178 was 
presented by Kellermann et al. (1998), as part of a Very Long 
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Baseline Array survey (VLBA 1 ) at 15 GHz. 

In this paper we present the first polarimetric VLBI obser- 
vations of 0735+178 at 1.3 cm and 7 mm, obtained with the 
VLBA. The images show an apparently quite twisted structure 
in the inner two milliarcseconds of the jet, in very good agree- 
ment with the total inten sity 2 cm VLBA image obtained by 
Kellermann et al. (1998). We analyze and interpret the mag- 



netic field struture and jet morphology, and discuss the possi- 
bility of a precessing jet in 0735+178. 

2. OBSERVATIONS 

The observations were performed on 1996 November 1 1 and 
December 22 using the VLBA at 1 .3 cm and 7 mm in snapshot 
mode. The data were recorded in 1-bit sampling VLBA format 
with 32 MHz bandwidth per circular polarization. The reduc- 
tion of the data was performed within the NRAO Astronomical 
Image Processing Syst em (A IPS) software in the usual manner 
(e.g., Leppanen et al. 1995). Delay differences between the 



right- and left-handed systems were estimated over a short scan 
of cross-polarized data of a strong calibrator (3C 454.3). The 
instrumental polarization was determined usi ng the feed solu- 
tion algorithm developed by Lepp anen et al. ( 1995[ ). We refer 
the readers to Gomez et al. (1998) for further details about the 



reduction and calibration of the data. 

Due to poor weather conditions in some of the antenna loca- 
tions during the second run, specially at Brewster, much of the 7 
mm data had to be deleted, leading to a significantly poorer total 
intensity image and the loss of all polarization at this epoch. 

3. RESULTS 

Figures [l] and || show the VLBA total and linearly polarized 
intensity images of 0735+178 at 1.3 cm and 7 mm, respectively. 
Tables [j] and g summarize the physical parameters obtained for 
0735+178 at 1.3 cm and 7 mm, respectively. Tabulated data 
corresponds to position in right ascension and declination rela- 
tive to the core component, total flux, polarized flux, magnetic 
field position angle, degree of polarization, and separation and 
structural position angle relative to the core component. Com- 
ponents in the total intensity maps were analyzed by model- 
fitting the uv data with circul ar Gau ssian components using the 
Difmap software (Shepherd |1997| ). Components are labeled 
from east to west using upper-case letters for the images at 1.3 
cm, and are marked in Figs. |IJ and [| Due to the potentially 
large uncertainties associated with the process of model-fitting, 
only components well above the noise level were considered. 
For example, there is some weak evidence for the existence of 
components d and cl in the 1 .3 cm maps, but the errors in the 
model-fitting are sufficiently large that we have decided not to 
include these in the final model. Component B is best fit at the 
second 7 mm epoch by two separate components labeled b2 and 
bl, the combined flux of which compares quite well with that 
of component b in the 7 mm image at the first epoch. To reduce 
the errors in the model-fitting we have fit the emission from the 
section of the jet beyond component B as a single component, 
labeled A. 

3.1. The Twisted Structure of the Inner Jet of 0735+178 



Baath et al. (|1991j) and Baath & Zhang ( |1991[ ) suggested the 
existence of bent structure in the jet of 0735+178 in order to ex- 
plain the apparent acceleration of superluminal components as 

1 The VLBA is an instrument of the National Radio Astronomy Observatory, 
agreement by Associated Universities, Inc. 



they moved further from the core. G94 found that components 
seem to move along linear trajectories, but to be ejected with 
different position angles, which could lead to apparent curved 
structure in the VLBI jet. 

Our observations directly confirm the existence of a very 
twisted structure in the inner region of the jet in 0735+178, 
in good agre emen t with the observations presented by Keller- 
mann et al. (1998). In our 1.3 cm maps we observe two sharp 
90° bends in the projected trajectory of the jet, one near the po- 
sition of component C and the other close to B. At 7 mm the 
bend near the region of component c2 (component C at 1.3 cm) 
is also quite evident, while the second bend can only be traced 
at the second epoch, since the emission beyond component b is 
resolved out in November 1996. 

Through a comparison with our maps, we can offer an ex- 
planation for the otherwise puzzling structure Zhang & Baath 
(1991 ) observed in their 1.3 cm VLBI maps. It is possible that 
they observed distinct components at different positions of a 
common curved path, consistent with that traced in our im ages. 
However in the 1.3 cm VLBI map by Baath et al. (1991) they 
observed in the inner milliarcsecond a straight jet extending 
northeast at a position angle of ~45°, in contrast with the struc- 
ture we observe in our images. Furthermore, G94 measured 
different ejection position angles for components in 0735+178. 

In order to allow an easier comparison with our images, we 
have marked in our November 1996 map of Fig. [j] the expected 
position for components K3, K4, K5, and K6 of G94 using 
their measured proper motions and assuming ballistic trajecto- 
ries from the core since their last observing epoch in 1992.44. 
For component K5 we have used the position reported for their 
1990.47 epoch, otherwise the extrapolation would place this 
component far too east of C, outside of the structure we have 
mapped. No proper motion data is available for component 
K6, but assuming a similar value to that observed for the other 
components, we have used 0.3 mas yr" 1 . Figure [l] shows that 
this extrapolation places the G94 components within the struc- 
ture found on our images, suggesting ballistic motion with sys- 
tematically changing component ejection directions as the most 
plausible explanation for the curved structure of 0735+178. 

However, the data is also consistent with motion of compo- 
nents following a common curved path. Component K3 was 
observed to change drastically its position angle relative to the 
core between 1987.41 and 1990.47, from 75° to 44° (see G94). 
Component K5 also experienced a change in its position an- 
gle between 1990.47 and 1992.44, with its velocity vector be- 
coming more aligned toward the direction of component C, be- 
tween the position of components c3 and c2. The positions 
of components G94 in the inner milliarcsecond are consistent 
with the twisted structure we have mapped, therefore we can- 
not rule out the possibility of a common curved funnel through 
which components flow. Note also that this curved funnel could 
change with time if it were produced by precession of the noz- 
zle in 0735+178. In the case of ballistic motion, the direction 
of ejection of new components should have changed progres- 
sively from 54°, to 35°, and then to 15° from 1984 to 1991.25 
to explain the different structural position angles observed for 
the birth of K4, K5, and K6. The direction of ejection would 
then have needed to change to 78° within about 4 yr, if we 
assume that component C was ejected around the beginning of 
1995, coincident with an outburst in total fl ux at 2 2 and 37 GHz 
measured at Metsahovi (Terasranta et al. 1998). Since then, 
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FIG. 1 . — VLBA images of 0735+178 at 1.3 cm at epochs 1996 November 1 1 (left) and 1996 December 22 (right). Total intensity is plotted as contour maps, 
while the linear gray scale shows the linearly polarized intensity. The superposed sticks give the direction of the magnetic field vector. Contour levels are -0.5, 0.5, 
1, 2, 4, 8, 16, 32, and 64% of the peaks of 335 (November 1996) and 322 (December 1996) mjy beam" 1 . For November 1996 the convolving beam is 0.74 X 0.38 
mas (FWHM), with position angle -1.5°, and the maximum in polarized intensity is 7.6 mJy beam -1 , with a noise level of 1.5 mjy beam -1 . For December 1996 the 
convolving beam is 0.73 X 0.36 mas (FWHM), with position angle -2.3° , and the maximum in polarized intensity is 12.9 mJy beam -1 , with a noise level of 2.6 mjy 
beam -1 . 




FIG. 2.— Same as Fig. |l|, but at 7 mm. Contour levels are -1, 1, 2, 4, 8, 16, 32, and 64% of the peaks of 307 (November 1996) and 362 (December 1996) 
mJy beam -1 . For November 1996 the convolving beam is 0.57 X 0.24 mas (FWHM), with position angle -20°, and the maximum in polarized intensity is 8.8 mJy 
beam -1 , with a noise level of 3.5 mjy beam -1 . For December 1996 the convolving beam is 0.93 X 0.43 mas (FWHM), with position angle -30.6° ; no polarization 
information was obtained. 



Table 1 

22 GHz models for 0735+178. 



a S I p Xm m r 

Component (mas) (mas) (mJy) (mJy) (°) (%) (mas) (°) 



1996.86 

Core 354 7 86 2 

C 0.81 0.19 361 12 -12 3 0.83 77 

B 0.7 1.3 87 9 31 10 1.5 29 

A 3.6 1.8 37 4.1 63 

1996.98 

Core 340 13 72 4 

C 0.82 0.18 304 9 -21 3 0.84 78 

B 0.6 1.2 94 3 36 3 1.3 29 

A 3.2 1.8 39 3.7 61 
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it would have had to remain at a similar orientation (at least 
on the plane of the sky) to give birth to components c3 and 
d, which according to two small flares in the Metsahovi data 
(Harri Terasranta, private communication) took place around 
1995.8 and 1996.7 (just before our first epoch), respectively. 
Note that the apparently curved jet of 0735+178 would ap- 
pear more pronounced if components are ejected with different 
speeds, as seems to be the case. 

Because of the short time range between our two epochs, the 
differences in the positions of the components lie within the 
errors in the model-fitting, therefore we can only rely on the 
observed magnetic field structure, and previous detections of 
components motions, to study the possibility of ballistic mo- 
tions for the components in 0735+178. 

3.2. Polarization 

The peak polarization at both epochs and wavelengths is lo- 
cated at the core, with a degree of polarization between 2 and 
4%, very close to the values measured by Gabuzda, Wardle, & 
Roberts (1989) and G94 at 6 and 3.6 cm. However, these values 



should be regarded as approximate, since maxima in the images 
of polarized intensity are not always precisely coincide nt with 
maxima in the total intensity (see also Gomez et al. 1998). 



Greater differences between the peaks in total and polarized in- 
tensity are found in the 7 mm map. Some polarized emission is 
also visible further from the core, especially at the first epoch 
at 1 .3 cm. In this map the jet between components C and B can 
be traced in polarization, where it has a rather uniform inten- 
sity. For the second epoch at 1.3 cm the core and component 
C are clearly detected in polarization, and there is also weaker 
polarized flux detected at the position of B. In the 7 mm po- 
larization map the core and component c2 are clearly visible. 
Another component appears close to the core, and we have ten- 
tatively identified this component as the polarized counterpart 
of d. Component c3 is also observed in polarization, while cl 
is just marginally detected. 

A very similar orientation of the magnetic field vector in the 
jet is observed at both epochs in the 1.3 cm maps and in the 
7 mm map at November 1996. The magnetic field in the core 
at our first epoch is oriented toward the east, in the direction 
of component C. This is cons istent with the values found by 
Gabuzda, Wardle, & Roberts (1989), although opacity effects 
due to the different wavelengths of observation may render this 
comparison meaningless. If the core is optically thick at both 
wavelengths, we should consider an extra rotation of 90° in the 
polarization angle when comparing to the remaining optically 
thin jet. At epoch November 1996, adding the flux of the core 
and component d at 7 mm, we obtain a very similar value (the 
same to within the uncertainty) as that measured for the core 
at 1.3 cm. The same is true for the second epoch. The core 
polarization angle in the 7 mm map shows a slightly different 
value than at 1.3 cm, while component d shows a magnetic field 
more aligned with the jet. These differences between the 1.3 cm 
and 7 mm maps are probably due to differences in resolution, 
although they may reflect a change in the opacity of the core. 
G94 found a rotation of the polarization angle in the core when 
observed at different epochs and frequencies. These changes in 
the polarization angle may be due to ejections of new compo- 
nents with polarization angles different than the core. Because 
of the lower resolution, these may appear in their maps as varia- 
tions of a single blended component -their core- similar to what 
we have experienced when comparing our 1.3 cm and 7 mm 
maps. 



The magnetic field structure in the jet between C and B is 
complex, and difficult to interpret unambiguously. One obvi- 
ous possibility is that the magnetic field follows the direction of 
the apparently twisted jet flow. In this case, the magnetic field 
throughout this part of the jet would be longitudinal. However, 
another interpretation is possible: that different features move 
from the core in different structural position angles, and that the 
magnetic field vectors are transverse to the flow direction in C 
and along the flow direction from the core in B. This alternative 
point of view is also consistent with the fact the magnetic field 
in component c3 also points back toward the core. In either of 
these interpretations the magnetic field in the jet north of C is 
longitudinal. 

4. THEORETICAL INTERPRETATIONS AND CONCLUSIONS 

Although 0735+178 represents the most abruptly bent jet ob- 
served at milliarcseconds scales to our knowledge, many jets in 
BL Lacs and quasars appear to be curved. While this is ampli- 
fied by projection effects, since relativistic jets in blazars must 
be pointing close to the line of sight to explain superluminal 
motion and related phenomena, it is not clear what causes the 
bends of at least 5° that must be present to account for the ob- 
served twists. 



One idea (Hardee 1987 ), whi ch has supporting observational 
evidence (Denn & Mutel 1998), is that the jet precesses, with 
the velocity vectors of the flow following the bends in the 
jet axis. Numerical three-dimensional magnetohydrodynamical 
simulations (Hardee, Clarke, & Rosen 1997) show how Kelvin- 



Helmholtz (K-H) helical instabilities can develop from an ini- 
tial induced precession at the jet inlet, resulting in a twisted jet 
with a helical geometry and magnetic field configuration. Such 
K-H instabilities can also arise from an external pressure gradi- 
ent not aligned with the initial direction of the jet flow. 

Another possibility for jet bending is ballistic precession, 
such that the flow velocity is directed radially away from the 
jet apex, in the manner of SS 433; this may be suggested by the 
trajectories of components K3, K4, K5, and K6 (G94). Alter- 
natively, the nozzle of the jet could change direction in a more 
erratic fashion than for precession. Bends can also occur in 
response to gradients in the pressure of the external medium 
that confines the jet. Another possibility, also suggested for 
0735+178 by G94, is that components only illuminate part of a 
broad jet funnel, giving the illusion of a bent trajectory. 

Each possibility for jet bending results in specific predic- 
tions. For example, in the case of erratic change in the nozzle 
direction, knots should separate from the core following bal- 
listic trajectories. For ballistic precession of the nozzle, the 
different sections of the jet should connect smoothly and the 
magnetic field direction should relate to the radius vector from 
the core, rather than follow the jet curvature. For precession 
in which the flow velocity vector changes direction with the jet 
axis, moving components should follow the bends, as should 
the magnetic field vector. In addition, the flux density and ap- 
parent velocity of the components should change as the angle 
between the vel ocity and the line of sight vary, as in 4C 39.25 
(Alberdi et al. |1993[ ). A precessing jet can be distinguished 
from one that changes because of external pressure gradients 
by the regular pattern of the bend, which should follow the ge- 
ometry of a helix in projection. 

With only 41 days between our two epochs, the expected 
motions of the components in the jet of 0735+178 lie within 
the errors in our model-fitting, hence we cannot yet conclude 
whether components on these scales follow ballistic trajecto- 
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Table 2 
43 GHz models for 0735+178. 



a 8 I p Xm tn r 8 

Component (mas) (mas) (mJy) (mJy) (°) (%) (mas) (°) 



1996.86 

Core 241 7 116 3 

d 0.10 0.02 98 2 84 2 0.10 80 

c3 0.46 0.11 26 4 29 15 0.47 76 

c2 0.87 0.09 189 10 -6 5 0.87 84 

cl 0.96 0.66 62 2 26 3 1.16 56 

b 0.7 1.5 48 1.7 24 

1996.98 

Core 370 

c 0.78 0.18 229 0.80 77 

b2 0.8 1.2 23 1.4 34 

bl 1.1 1.7 23 2.0 33 

a 3.4 1.9 13 3.9 61 



ries or the direction of the jet axis. A ballistic extrapolation of 
the position of the components observed in G94 shows that bal- 
listic motions could explain the observed twisted structure for 
the jet. The information provided by the direction of the mag- 
netic field, which follows the bend near component B, suggests 
that the jet in 0735+178 precesses such that the flow velocity is 
parallel to the jet axis. Indeed, the curved jet axis is consistent 
with the projection of a helix, which would have to experience 
an increase in the helical wavelength beyond component B to 
account for the rather straight jet observed at larger scales (see 
G94). However, the interpretation of the magnetic field struc- 
ture in B and in the jet between B and C is not entirely clear. It 
may be that the magnetic field in this part of the jet is associated 
with radial flow from the core along different position angles, 
and is longitudinal -directed back toward the core, suggestive 



of ballistic motion. 

Further high resolution polarimetric VLBI observations of 
0735+178 are necessary to measure proper motions and deter- 
mine whether the apparent twisted structure in this source rep- 
resents a common curved path for the jet flow or a superposi- 
tion of independently ejected ballistically moving components. 
Such observations would also make it possible to search for sys- 
tematic changes in the jet emission of the sort expected if the 
nozzle of the jet of 0735+178 precesses. 
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